Combined effects of thermal radiation and Hall 
current on MHD free-convective flow and mass 
transfer over a stretching sheet with variable 

viscosity 

G. C. Shit* and R. Haldar 

Department of Mathematics 
Jadavpur University, Kolkata - 700032, India 

Abstract 

An analysis has been investigated for the effects of thermal radiation and Hall current 
on magnetohydrodynamic free-convective flow and mass transfer over a stretching sheet 
with variable viscosity in the presence of heat generation/absorption. The fluid viscosity 
is assumed to vary as an inverse linear function of temperature. The boundary-layer equa- 
tions governing the flow problem under consideration have been reduced to a system of 
non-linear ordinary differential equations by employing a similarity transformation. Us- 
ing the finite difference scheme, numerical solutions to the transform ordinary differential 
equations have been solved and the results that obtained are presented graphically. With 
an aim to test the accuracy, the numerical results have been compared with the existing 
scientific literature and found excellent agreement. 

Keywords: Thermal radiation; Variable viscosity; MHD; Hall current; Heat and Mass 
transfer 

1 Introduction 

In the recent past there has been a growing interest in boundary-layer flow on a continuous 
moving surface in the presence of magnetic field with or without considering the effect of Hall 
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current. But, of course these are very significant type of flow towards the several engineering 
applications, such as in polymer processing, electro-chemistry, MHD power generators, flight 
magnetohydrodynamics as well as in the field of planetary magnetosphere, aeronautics and 
chemical engineering. 

Sakiadis [1] first explored the study of boundary-layer flow on a continuous moving surface 
and Crane [2] extended this problem to a stretching sheet whose surface velocity varies linearly 
with the distance x from a fixed origin. During the past decades several investigators [3-10] 
have considered the boundary- layer flow problems under different physical situations. Gupta 
and Gupta [11] examined the heat and mass transfer over a stretching sheet subject to suction 
or blowing. The influence of a uniform magnetic field on the flow of an electrically conduct- 
ing fluid past a stretching sheet was investigated by Pavlov [12], Andersson [13], Andersson et 
al.[14], Gupta and Chakrabarty [15], Char [16], Watanabe and Pop [17] and Elbashbeshy [18]. 
The chemical reaction on free-convective flow and mass transfer of a viscous, incompressible and 
electrically conducting fluid over a stretching sheet was investigated by Afify [19] in the presence 
of a uniform transfer magnetic field. In all these investigations the electrical conductivity of the 
fluid was assumed to be uniform and low magnetic field intensity. However, in an ionized fluid 
where the density is low and thereby magnetic field intensity is very strong, the conductivity 
normal to the magnetic field is reduced due to the spiraling of electrons and ions about the 
magnetic lines of force before collisions take place and a current induced in a direction normal 
to both the electric and magnetic fields. This phenomena is known as Hall effect. Watanabe 
and Pop [20] investigated the magnetohydrodynamic boundary-layer flow over a continuously 
moving semi-infinite flat plate by taking into account the Hall currents. Aboeldahab [21] and 
Aboeldahab and Elbarbary [22] studied the Hall current effects on MHD free-convective flow 
past a vertical plate with mass transfer. Shit [23] investigated the Hall effects on MHD free- 
convective flow and mass transfer over a stretching sheet in the presence of chemical reaction. 
Fakhar et al. [24] studied the Hall effects on the unsteady magnetohydrodynamic flow of a 
third grade fluid without considering the heat and mass transfer phenomena. The effect of Hall 
currents on the steady MHD flow of Berger's fluid between two parallel electrically insulating 
infinite planes was carried out by Rana et al. [25]. 



2 



Recently, a new idea is added to the study of boundary-layer fluid flow and heat transfer is 
the consideration of the effect of thermal radiation and temperature dependent viscosity. Many 
processes in engineering applications occur at high temperatures and the radiate heat transfer 
becomes very important for the design of the pertinent equipment. In view of this, Rap its and 
Perdikis [26] and Rapits [27] studied respectively the flow of a visco-elastic fluid and micropolar 
fluid past a stretching sheet in the presence of thermal radiation. Mukhopadhaya et al. [28] in- 
vestigated the problem of M HD boundary-layer flow over a heated stretching sheet with variable 
viscosity. The radiation effect on boundary-layer flow with or without applying magnetic field 
under different situations were studied by Shateyi [29], Mahmoud [30], Pal and Talukdar [31] 
and Pal and Chatterjee [32]. However, Salem [33] investigated the effect of variable viscosity 
on MHD viscoelastic fluid flow and heat transfer over a stretching sheet without considering 
thermal radiation effect. Shit and Haldar [34] carried out the study of the effect of thermal 
radiation on MHD viscoelastic fluid flow past a stretching surface with variable viscosity. But, 
no attempt is available in the existing scientific literatures for the consideration of the combined 
effect of thermal radiation and Hall current on the study of MHD boundary-layer flow. Thus 
the present study fills the gap in this directions. 

Again, combined heat and mass transfer problems with chemical reaction are of increasing 
importance in many processes, like drying, evaporation at the surface of a water body, energy 
transfer in a wet cooling tower etc. In this context, Muthucumarswamy and Ganesan [35] stud- 
ied the effect of the chemical reaction and injection as well as flow characteristics in an unsteady 
upward motion of an isothermal plate. Chamakha [36] carried out the MHD flow of uniformly 
stretching vertical permeable surface in the presence of heat generation / absorption along with 
the chemical reaction. Very recently, Mohamed and Abo-Dahab [37] have investigated the in- 
fluence of chemical reaction and thermal radiation on hydromagnetic free-convective heat and 
mass transfer for a micropolar fluid via a porous medium bounded by an infinite vertical porous 
plate in the presence of heat generation. Seddeek et al. [38, 39] analyzed the effects of chemical 
reaction, radiation and variable viscosity on hydromagnetic mixed convection heat and mass 
transfer. In all these studies ignores the effect of the consideration of Hall current. 

Now we propose to study the combined effects of thermal radiation and Hall current on 
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the hydromagnetic free-convective flow and mass transfer over a stretching surface with variable 
viscosity in the presence of heat generation/absorption. The present problem pertains to situa- 
tion in which the n-th order chemical reaction takes place. Thus the study is also applicable to 
the elongation of the bubbles and in bioengineering where the flexible surfaces of the biological 
cells and membranes in living systems are typically surrounded with fluids which are electrically 
conducting and being stretched constantly. 

2 Mathematical Formulation 

We consider the steady free-convective flow and mass transfer of an incompressible, viscous and 
electrically conducting fluid past a stretching sheet and the sheet is stretched with a velocity 
proportional to the distance from a fixed origin O (cf. Fig. 1). A uniform strong magnetic field 
of strength B is imposed along the y-axis and the effect of Hall currents is taken into account. 
Taking Hall effects into account the generalized Ohm's law [40] may be put in the form : 



where represents the velocity vector, is the intensity vector of the electric field, is the 
magnetic induction vector, v e the magnetic permeability, the electric current density vector, 
m = ^ is the Hall current parameter, a the electrical conductivity, e the charge of the electron 
and n e is the number density of the electron. A very interesting is that effect of Hall current gives 
rise to a force in the z-direction which in turn produces a cross flow velocity in this direction 
and thus the flow becomes three-dimensional. 

The temperature and the species concentration are maintained at a prescribed constant values 
T w , C w at the sheet and and are the fixed values far away from the sheet. Since the 
concentration of diffusing species is very small in comparison to other chemical species, the Soret 
and Dufour effects are neglected. 

Following Lai and Kulacki [41], the fluid viscosity \x is assumed to vary as a reciprocal of a linear 
function of temperature given by 




4 



^ = ^-[1 + 7^-^)] (1) 

or 

— = a{T- T r ) (2) 

A^oo 

where a = — and T r = — -. 

Moo 7 

In the above equation both a and T r are constants, and their values depend on the thermal 
property of the fluid, i.e., 7. In general a > represent for liquids, whereas for gases a < 0. 
By assuming Rosseland approximation for radiation, the radiative heat flux q r is given by 

4a* <9T 4 

where a* is the Stefan-Boltzman constant and K the mean absorption coefficient. We assume 
that the temperature differences within the flow are sufficiently small such that T 4 may be 
expressed as a linear function of the temperature as shown in Chamakha [36]. Expanding T 4 in 
a Taylor series about and neglecting higher order terms, we obtain 

T 4 - ATlT - 3T 4 (4) 

Substituting T 4 from (4) in (3) and differentiating the resulting equation with respect to y, we 
obtain as 

dy 3K dy 2 1 } 

Owing to the above mentioned assumptions, the boundary layer free-convection flow with mass 
transfer and generalized Ohm's law is governed by the following system of equations : 



du Ou \ d ( du \ 
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( dw dw\ d ( dw\ oBl . , . . 



dc ac „a 2 c , ,„ r , ,„ 

where (u, v, w) are the velocity components along the (x, y, z) directions respectively, \i is the co- 
efficient of viscosity, g the acceleration due to gravity, j3 t the coefficient of thermal expansion, j3 c 
the coefficient of expansion with concentration, T and C are the temperature and concentration 
respectively, D the thermal molecular diffusivity, k is the constant measures the rate of reaction, 
o is the electrical conductivity, c p is the specific heat at constant pressure, k is the thermal con- 
ductivity, Too and poo are the free stream temperature and density and n be the order of reaction. 

The boundary conditions to the present problem can be written as 

u = bx, v = w = 0, T = T W , C = C W at y = (11) 

u — y 0, w ->■ 0, T — >■ T^, C ->■ Coo as y ->■ oo (12) 

where 6(> 0) being stretching rate of the sheet. The boundary conditions on velocity in (11) are 
the no-slip condition at the surface y — 0, while the boundary conditions on velocity at y — > 
follow from the fact that there is no flow far way from the stretching surface. 

To examine the flow regime adjacent to the sheet, the following transformations are invoked 
u = bxf'(rj); v = -Vbvf(rj); w = bxg(rj); rj = J-y; 6{rj) = ^ _" Z °° ; <f>{r}) = ^ (13) 

V ^ -L w 1 oo ^oo 

where / is a dimensionless stream function, rj is the similarity space variable, 9 and (f> are the 
dimensionless temperature and concentration respectively. Clearly, the continuity equation (6) 
is satisfied by u and v defined in (13) on substitution which into equations (7) - (10) gives 
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^) (/ ;_ /r)+r _(_i_) r _ (^*) 



(37Vr + 4)0" + 3NrPrf9' + 3NrPr\9 = 0, (16) 

^ + ScW-70 B )=O, (17) 
and the transformed boundary conditions are given by 

/'(?/) = 1, /(?/)= 0, 0(77) = 0, 0(t/) = 1, 0(r/) = l at 77 = 0, (18) 

f(77) = 0, (7(77) = 0, 0(77) = 0, ^(r/) = at 77 00, (19) 

where primes denotes differentiation with respect to 77 only and the dimensionless parameters 
appearing in the equations (14)-(17) are respectively r = = — ^ T J_ Toa ) is known as 

the viscosity parameter, M = ^| the magnetic parameter, P r = poo< ^ pV the Prandtl number, 
m = is the Hall current parameter^ = ^ (C w — Coo)" -1 the non-dimensional chemical 
reaction parameter, G r = pocg/3t g^~ Too) the Grashof number, G c = p^MCu,-c xj ) thg mo difi ec i 
Grashof number, Nr = J^ a , the thermal radiation parameter, Sc — the Schmidt number 
and A = p Q c b is defined as the heat generation/absorption parameter. 



The important characteristics of the present investigation are the local skin-friction coefficient 
Cf, the local Nusselt number Nu and the local Sherwood number Sh defined by 

C f = — ^ = r(0), where r w = y. (^\ = ^x\T-f\<S), (20) 



Nu= qw = -6'(0), where q w = -k(?f) — —k\ —(T w — T OO )0'(O), (21) 



V 



Sh= ™ w = - (f) '(Q), where m w = -D(^-\ = —D\i^-(C W — C OO )0'(O), (22) 



y=o 



If we consider M = m = and Nr = Sc = G r = G c = 0, the present flow problem becomes 
hydrodynamics boundary-layer flow past a stretching sheet whose analytical solution put for- 
warded by Crane [2] as follows : 



/fa) = l- e -> i.e, f( V ) = e-v (23) 

An attempt has been made to validate our results for the axial velocity f'(rj), we compared our 
results with this analytical solution and have found excellent agreement. 



3 Numerical Results and Discussion 

The system of coupled and non-linear ordinary differential equations (14)-(17) along with the 
boundary conditions (18) and (19) have been solved numerically by employing a finite difference 
scheme. We used Newton's linearization method (cf. Cebeci and Couteix [42]) to linearize the 
discretized equations. The essential features of this technique is that it is based on a finite dif- 
ference scheme, which has better stability, simple, accurate and more efficient. Finite difference 
technique leads to a system which is tri-diagonal and therefore speedy convergence as well as 
economical memory space to store the coefficients. The computational work has been carried 
out by taking Si] = 0.0125 and further reduction in 5r] does not bring about any significant 
change. In the present study, the numerical values to the physical parameters have been chosen 
so that M, m, Nr, 9 r , Pr, Gr, Gc, Sc,n,j and A are varied over a range, which are listed in the 
figure legends. Fig. 2 shows that our numerical results are complete agreement with those of 
Crane [2]. 
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Figs. 3-10 illustrate the variation of axial velocity for different values of the dimensionless 
parameters that involved in the present study. Fig. 3 shows that the axial velocity decreases 
with the increase of the magnetic parameter M, whereas from Fig. 4 it indicates that the axial 
velocity increases with the increase of Hall parameters m. This is due to the fact that as M 
increases, the Lorentz force which opposes the flow and leads to deceleration of the fluid motion. 
By contrast, the cross flow velocity component g(rj) induced due to Hall effects and shows a 
anomalous behaviour in f'{rj) with the variation of M. It has been shown in Figs 5 and 6 that 
the axial velocity decreases with the increase of the Prandtl number Pr as well as the thermal 
radiation parameter Nr. This is due to fact that there would be a decrease of boundary- layer 
thickness in the presence of thermal radiation. Fig. 7 depicts that the axial velocity f'(rj) in- 
creases with the decreasing of the viscosity parameter 9 r . This observation leads to an increase 
of the thermal boundary- layer thickness. The effects of heat generation parameter (A > 0) and 
absorption parameter (A < 0) on the axial velocity displayed in Fig. 8. This figure shows that 
the axial velocity decreases as the parameter A increases. The variation of Schmidt number Sc 
and the chemical reaction parameter 7 on the axial velocity profile f'(rj) shown in Figs. 9 and 
10 respectively. It is obvious that the increased values of Sc and 7 tend to decreasing of the 
velocity profiles across the boundary-layer. 

Figs. 11-19 give the distribution of the z- component of velocity, which is induced due to 
the presence of Hall effects. All these figures show that for any particular values of the physical 
parameters g{rj) reaches a maximum value at a certain height rj above the sheet and beyond 
which g{rj) decreases gradually in asymptotic nature. It is noticed from Fig. 11 is that in the 
absence of magnetic parameter M(= 0) cross flow velocity vanishes. This is due to fact that, 
when there is no applied magnetic field, the cross flow velocity would not arise. The variation of 
Hall current parameter m on the cross-flow velocity g(rj) shown in Fig. 12. An interesting result 
observed from this figure that the cross-flow velocity gradually increases with the increase of 
m < 2 and the velocity decreases for m > 2. The values of m beyond which the flow behaviour 
changes is considerable depend upon the choice of the magnetic parameter M. Thus we conclude 
that after certain magnetic field strength the flow behaviour is significantly affected. Figs. 13 
and 14 indicate that the cross- flow velocity g(rj) decreases with increasing the Prandtl number 
Pr and the thermal radiation parameter Nr, while from Fig. 15 that the trend is reversed after 
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a certain height above the sheet. This is lies in the fact that the increase of Prandtl number Pr 
and the thermal radiation Nr give rise to decrease of the momentum boundary-layer thickness. 
We observed from Figs. 16, 17 and 19 that the cross-flow velocity decreases with the increase 
of the heat generation/absorption parameter A, Schmidt number Sc and the chemical reaction 
parameter 7. However, the cross-flow velocity increases as the order of the chemical reaction n 
increases. 

The distribution of dimensionless temperature 6(rj) along the height of the stretching sheet 
for different values of the dimensionless parameters involved in the present study displayed 
through Figs. 20 - 25. Fig. 20 shows that by the application of an external magnetic enhances 
the temperature of the fluid, while the effect of the Hall current parameter m has an reducing 
effect on the dimensionless temperature 6{rf) shown in Fig. 21. It may note that the effect of 
Hall current parameter m opposes the effect of magnetic field on the temperature distribution. 
Fig. 22 presents the variation of Prandtl number Pr on the temperature 9(rf) . The results pre- 
sented in Fig. 22 shows that the dimensionless temperature decreases as the Prandatl number 
Pr increases. This is lies in the fact that smaller values of Pr are equivalent to increasing of 
thermal conductivities and therefore heat is able diffuse away from the stretching sheet. Fig. 
23 represents the temperature profiles for various values of the thermal radiation parameter 
Nr in the boundary- layer. Increasing the thermal radiation parameter Nr produces a decrease 
in the temperature of the fluid. This is because of the fact that the thermal boundary-layer 
thickness decreases with increasing the thermal radiation parameter. Fig. 24 gives the variation 
of the viscosity parameter 9 r on the temperature profiles and which shows that no significant 
change occur with the increase of the values of 6 r . It is evident that the energy equation (16) 
is uncoupled from the viscosity parameter 6 r . For this reason no figures are presented herein 
for the variation of Sc, n and 7. The effect of heat generation/absorption parameter A on the 
dimensionless temperature 9{rj) is shown in Fig. 25. It is clear that an increase in the heat 
generation/absorption parameter A leads to a decrease of #(77). Thus the effect of internal heat 
generation is to decrease the rate of energy transport to the fluid, thereby decreasing the tem- 
perature of the fluid. 

The effect of the imposition of various parameters on concentration profiles are shown in 
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Figs. 26-31. It is observed from Figs. 26 and 27 that under the action of a strong mag- 
netic field, the concentration species has an enhancing effect, whereas it has the reducing effect 
on the Hall current parameter m. The reason behind this is that in equations (14) and (15), 
the parameters M and m are connected by the relations of the form and ■ As the 
variation of concentration profiles for different values of the parameters Nr and 6 r agrees with 
the temperature profiles, we do not have presented those results. However, Fig. 28 is the vari- 
ation of heat generation/absorption parameter A on the species concentration, which opposes 
the variation of the temperature distribution. Fig. 29 depicts the variation of Schmidt number 
Sc on the species concentration <f)(rf). It is observed that the species concentration decreases 
with the increase of the Schmidt number Sc. Physically, which shows that the increase of Sc 
causes decrease of molecular diffusion D. The influence of chemical reaction parameter 7 on 
the species concentration profiles for generative chemical reaction is shown in Fig. 30. It is 
noticed that the increasing values of the chemical reaction parameter 7 is to decelerate the 
concentration distribution in the boundary- layer. This is due to fact that destructive (7 > 0) 
chemical reaction decreases the shortest boundary-layer thickness and thereby increasing the 
species concentration. It is observed from Fig. 31 that the increasing effect of the order of 
chemical reaction parameter n is to enhances the mass transfer in the boundary- layer. 

The problems of engineering interest for the present study are the local skin-friction coef- 
ficient Cf , the local Nusselt number Nu and the local sherwood number Sh which indicates 
physically the wall shear stress, the rate of heat transfer at the sheet and the rate of mass trans- 
fer respectively. The expressions of these physical quantities have been presented in equations 
(20)- (22). Tables 1-3 exhibit the numerical values to the local skin-friction coefficient /"(0), 
the local Nusselt number — 0'(O) and the local sherwood number — 0'(O) respectively. It has 
been observed that for any particular values of the parameters Pr, Nr, m, 9 r , 7 , A and Sc, 
the local skin-friction coefficient, the local Nusselt number and the local Shearwood number 
decreases with the increase of the magnetic parameter M. Table -1 shows that for fixed value of 
M, the local skin- friction Cf decreases with increasing the values of the Prandatl number Pr, 
thermal radiation parameter Nr, viscosity parameter r , the rate of chemical reaction 7, heat 
generation/absorption parameter A and the Schmidt number Sc. Physically, we meant that for 
the increasing of the thermal radiation parameter Nr leads to decreasing in the boundary-layer 
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thickness. However, the skin-friction coefficient increases as the Hall parameter m increases. 
It is worthwhile to observed from Table-2 that the rate of heat transfer at the sheet increases 
with the increase of the parameters Pr, Nr, m and 7. The increase of the thermal radiation 
parameter Nr has an enhancing effect on the thermal boundary-layer thickness. It is also no- 
ticed that the increasing of the parameters r , 7 and Sc is to reduce the heat transfer rate and 
which shows that the rate of change of heat transfer is insignificant. The rate of mass transfer 
increases significantly when the Hall parameter m, the chemical reaction parameter 7 and the 
Schmidt number Sc increases. Moreover, the rate of mass transfer decreases insignificantly with 
the increase of the parameters Pr, Nr and A. 

4 Conclusions 

In the present investigation, we dealt with the combined effects of thermal radiation and Hall 
current on the boundary-layer fluid flow, heat and mass transfer with temperature dependent 
viscosity. The highly non-linear coupled system of partial differential equations characterizing 
the flow, heat and mass transfer has been reduced to a coupled system of non-linear ordinary 
differential equations by applying a suitable similarity transformations. The resulting system 
solved numerically by using the finite difference scheme along with the Newton's linearzation 
technique. The obtained numerical results have been presented through the figures and in tab- 
ular form to illustrate the details of the flow behaviour, heat and mass transfer phenomena and 
their dependence on the physical parameters that involved in the present investigation. From 
our computed numerical results we observed that the magnetic field and Hall current produce 
opposite effects on the velocity distribution and heat transfer as well as on the concentration 
distribution. For a fixed value of M, the skin-friction increases with an increase in m and similar 
is the observation for heat and mass transfer rate. Sufficiently strong heat generation parameter 
may alter the temperature gradient. Temperature decreases and the concentration increases 
with increasing values of Prandtl number Pr. But a reversal trend is observed when the val- 
ues of the thermal radiation parameter increases. The species concentration decreases with an 
increase in the values of the Schmidt number Sc and chemical reaction parameter A whereas 
opposite trend is observed in the case of the order of the chemical reaction. It is hopped that 



12 



the results obtained will serve as a scientific tool for understanding more complex flow problems 
and provide more useful information for engineering applications. 
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Fig. 1 Physical sketch of the problem 
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Fig. 6 Variation of firf) with r\ for different values of Nr 
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Fig. 10 Variation of f'(r)) with r] for different values of 7 
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12 Variation of g(rj) for different values of the Hall parameter m. 



22 



0.1 



gen) 




Fig. 13 Variation of g{rj) with rj for different values of Pr 
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Fig. 15 Variation of g(rj) with 77 for different values of 9 r 
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Fig. 16 Variation of g(rj) with 77 for different values of A 
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Fig. 17 Variation of g{rj) with rj for different values of Sc 
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Fig. 18 Variation of g(rj) with i] for different values of n 
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Fig. 21 Distribution of dimensionless temperature 6{rf) with rj for different values of m 
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Fig. 22 Distribution of dimensionless temperature 9{rf) with r) for different values of Pr 
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Fig. 25 Distribution of dimensionless temperature 6{rj) with rj for different values of A 
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Fig. 26 Influence of concentration species <p(r]) with r\ for different values of M 
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27 Influence of concentration species <j){rj) with rj for different values of m 
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Fig. 30 Influence of concentration species 4>{rf) with r\ for different values of 7 
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Fig. 31 Influence of concentration species (f>(rj) with rj for different values of n 
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Table-1. Numerical values of the local skin-friction coefficient Cf = f"(0) 



Pr Nr m 9 r 7 A Sc 


M = 0.0 


M = 2.0 


M = 4.0 


0.71 1.0 0.5 -2.0 0.5 0.5 0.5 
1.0 
3.0 


-0.52597 
-0.55356 
-0.66289 


-1.45996 
-1.47541 
-1.55043 


-2.10233 
- 2.11350 
-2.17159 


0.71 0.5 0.5 -2.0 0.5 0.5 0.5 
1.0 
5.0 


-0.49627 
-0.52596 
-0.57792 


-1.44462 
-1.45996 
-1.49009 


-2.09153 

- 2.10233 

- 2.12436 


0.71 1.0 1.0 -2.0 0.5 0.5 0.5 
3.0 
5.0 


-0.52596 

- 0.52596 

- 0.52596 


-1.22516 
- 0.73786 
-0.61575 


-1.76392 
-1.00484 
- 0.75757 


0.71 1.0 0.5 -1.0 0.5 0.5 0.5 
-2.0 
-5.0 


-0.56098 
-0.52596 
-0.49443 


-1.66484 
- 1.45996 
-1.31365 


-2.41425 

- 2.10233 

- 1.88306 


0.71 1.0 0.5 -2.0 0.1 0.5 0.5 

0.5 
1.0 


-0.49746 
-0.52596 
- 0.55088 


-1.44034 
-1.45996 
-1.47691 


-2.08909 
-2.10233 
-2.11402 


0.71 1.0 0.5 -2.0 0.5 -0.5 0.5 

0.0 
0.5 


-0.42172 
-0.48498 
-0.52596 


-1.38284 
-1.43029 
- 1.45996 


-2.04079 
-.07862 
-2.10233 


0.71 1.0 0.5 -2.0 0.5 -0.5 0.5 

1.0 
1.5 


-0.52596 
-0.57899 
- 0.61068 


-1.45996 
-1.48921 
- 1.50989 


-2.10233 
-2.12044 
-2.13432 
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Table-2. Numerical values of the local Nusselt number Nu = —0'(0) 



Pr Nr m 9 r 7 A Sc 


M = 0.0 


M = 2.0 


M = 4.0 


0.71 1.0 0.5 -2.0 0.5 0.5 0.5 
1.0 
3.0 


0.51857 
0.61651 
1.09706 


0.48436 
0.57475 
1.03107 


0.46751 
0.55368 
0.99413 


0.71 0.5 0.5 -2.0 0.5 0.5 0.5 
1.0 
5.0 


0.41891 
0.51857 
0.70871 


0.39371 
0.48436 
0.66088 


0.38159 
0.46751 
0.63618 


0.71 1.0 1.0 -2.0 0.5 0.5 0.5 
3.0 
5.0 


0.51857 
0.51857 
0.51857 


0.48973 
0.50769 
0.51382 


0.47213 
0.49279 
0.50509 


0.71 1.0 0.5 -1.0 0.5 0.5 0.5 
-2.0 
-5.0 


0.51673 
0.51857 
0.52016 


0.48081 
0.48436 
0.48729 


0.46364 
0.46751 
0.47070 


0.71 1.0 0.5 -2.0 0.1 0.5 0.5 

0.5 
1.0 


0.52097 
0.51857 
0.51662 


0.48641 
0.48436 
0.48279 


0.46901 
0.46751 
0.46639 


0.71 1.0 0.5 -2.0 0.5 -0.5 0.5 

0.0 
0.5 


0.12334 
0.35575 
0.57857 


0.02936 
0.30484 
0.48436 


-0.02061 
0.27964 
0.46751 


0.71 1.0 0.5 -2.0 0.5 -0.5 0.5 

1.0 
1.5 


0.51857 
0.51425 
0.51223 


0.48436 
0.48156 
0.48004 


0.46751 
0.46573 
0.46469 
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Table-3. Numerical values of the local Sherwood number Sh = —<f>'(0) 



Pr Nr m 9 r 7 A Sc 


M = 0.0 


M = 2.0 


M = 4.0 


0.71 1.0 0.5 -2.0 0.5 0.5 0.5 
1.0 
3.0 


0.67199 
0.66905 
0.66032 


0.62481 
0.62322 
0.61772 


0.60072 
0.59979 
0.59634 


0.71 0.5 0.5 -2.0 0.5 0.5 0.5 
1.0 
5.0 


0.67544 
0.67199 
0.66670 


0.62659 
0.62481 
0.62188 


0.60175 
0.60072 
0.59898 


0.71 1.0 1.0 -2.0 0.5 0.5 0.5 
3.0 
5.0 


0.67199 
0.67199 
0.67199 


0.63256 
0.65744 
0.66566 


0.60770 
0.63733 
0.65406 


0.71 1.0 0.5 -1.0 0.5 0.5 0.5 
-2.0 
-5.0 


0.66593 
0.67199 
0.67411 


0.61956 
0.62481 
0.62911 


0.59486 
0.60072 
0.60553 


0.71 1.0 0.5 -2.0 0.1 0.5 0.5 

0.5 
1.0 


0.50642 
0.67199 
0.83627 


0.43681 
0.62481 
0.80234 


0.40050 
0.60072 
0.78493 


0.71 1.0 0.5 -2.0 0.5 -0.5 0.5 

0.0 
0.5 


0.68296 
0.67616 
0.67199 


0.63338 
0.62786 
0.62481 


0.60670 
0.60279 
0.60072 


0.71 1.0 0.5 -2.0 0.5 -0.5 0.5 

1.0 
1.5 


0.67199 
0.96662 
1.19569 


0.62481 
0.90329 
1.12354 


0.60072 
0.86889 
1.08281 
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